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Abstract. We investigate the equilibration and thermalization properties of quantum 
systems interacting virith a finite dimensional environment. By exploiting the concept 
of time averaged states, we introduce a completely positive map which allows to 
describe in a quantitative way the dependence of the equilibrium state on the initial 
condition. Our results show that the thermalization of quantum systems is favored if the 
dynamics induces small system-environment correlations, as well as small changes in the 
environment, as measured by the trace distance. 



1. Introduction 

The mechanisms behind thermahzation have recently attracted a renewed interest and 
initiated the development of novel statistical formulations of equilibration in the realm of 
quantum mechanics [H [2 El HI El E H El E]- In all these descriptions the total system 
under investigation is associated with a finite dimensional Hilbert space and thus the 
asymptotic limit of the dynamics does not exist and the system returns, with possibly 
very long recurrence time, arbitrarily close to its initial state infinitely many times pUl lllj . 
Relaxation to equilibrium in the usual sense is thus impossible. Nevertheless, one can 
introduce an extended notion of equilibration if the system tends towards some state, which 
can be identified as equilibrium state of the dynamics, and stays close to it most of the time. 
There will still be some fluctuations around the equilibrium state, but extremely small or 
rare. 

In this work we consider the situation in which a closed quantum system can be 
decomposed into two parts, an open system S and a bath B, and investigate the equilibration 
properties of the subsystem S. We use the extended notion of equilibration, i.e, we will say 
that the open system equilibrates if its time evolved state (also called reduced state) ps{t) 
approaches some equilibrium state and spends most of the time close to it. In the same spirit, 
one can introduce the notion of thermalization of the open system by means of additional 
conditions on its equilibrium state [5] . Namely, one requires that the latter does not depend 
on the initial total state, beside a possible dependence on macroscopic parameters, such 
as temperature, characterizing the initial state of the bath. In this case one says that the 
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open system thermalizes if, in addition, the equilibrium state takes the form of a Gibbs 
state. The capabiUty of an open system to thermaUze ultimately traces back to specific 
properties of the total Hamiltonian, which fixes the evolution of the total closed system and, 
in particular, characterizes the interaction between the open system and the bath. Indeed, 
if the open system and the bath do not interact no thermalization is expected. Moreover, 
if there are conserved quantities for the open system its equilibrium state will unavoidably 
depend on the initial reduced state. More generally, it has been shown [8] that the lack of a 
sufficient amount of entanglement in the energy eigenbasis is a basic reason for the absence 
of thermalization. 

Here, we want to discuss the thermalization of open quantum systems within the above- 
mentioned framework, with the aim of clarifying the role of general dynamical mechanisms 
that can induce or prevent it. Apart from the obvious situation of conserved quantities for 
the open system, which features of the dynamics imply a dependence of the equilibrium state 
on the initial reduced state? Thermalization requires a sufficient amount of entanglement 
in the energy eigenbasis, but what is the role played by the interaction induced correlations 
between the open system and the bath? We investigate how the equilibrium state is 
modified when one resets the correlations between the system and the bath, as well as 
the environmental state, to their initial value. We find that small system-environment 
correlations, together with small changes in the environmental state, generally lead to an 
equilibrium state that hardly depends on the initial state of the open system. The system- 
environment correlations built by the dynamics thus play a role which is in a sense opposite 
to that of the entanglement in the energy eigenbasis. The interaction induced correlations 
can prevent the thermalization: strong system-environment correlations and changes in 
the environmental state allow to distinguish between equilibrium states corresponding to 
different initial states. The information about the initial state of the reduced system, 
transferred to the environment through the establishment of correlations, does infiuence 
the equilibrium state. In this sense it can be considered as trapped in the equilibrium state. 
For this reason, we will refer to this mechanism preventing thermalization as information 
trapping. The quantitative characterization of information trapping will be given in terms 
of trace distance, which measures the distinguishability between quantum states [12] and 
has already been used to detect through its variation the information flow between system 
and environment [HI \W\ , as also discussed later on. 

2. General framework 

Consider a finite dimensional Hilbert space, which can be decomposed as "H = 'Hs'Si'Hb, with 
Hs and Hb Hilbert spaces associated with the open system S and the bath B, respectively. 
A crucial assumption here is that the total Hilbert space H has finite dimension, which will 
be denoted by d, while we denote by ds and ds the dimensions of the open system and the 
bath. The evolution of the total system is governed by a one-parameter group of unitary 
operators U (t) , which is fixed by the total Hamiltonian 

H = Hs + Hb+Hsb. (1) 

where Hg and Hb are the free Hamiltonian of the system and the bath, and HgB is the 
interaction term. Moreover, assume a product initial total state, so that, for a fixed initial 
state of the bath, there is a well-defined reduced dynamics [T3], i.e. there is a family of 
completely positive and trace preserving maps A{t) on the set S{Hs) of statistical operators 
on Us such that the reduced state psit) at time t is given by ps{t) = A{t)ps, where 
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PS = Ps(0): 

A{t) : S{Hs) SiHs) 

PS ^ Psit)=A{t)ps. (2) 

In [3j it has been shown that an open system equihbrates under very general assumptions 
if the effective dimension of the bath, i.e. the dimension of the subspace of involved 
into the dynamics, is much larger than the open-system dimension ^5. In this case, for any 
initial state ps, the corresponding time evolved state ps{t) will be most of the time close to 
the time averaged state ps, which is defined as 

_ 1 /•* 

PS = hm - / drpsir) (3) 

t-i-OO t Jq 

and represents the equilibrium state of the reduced dynamics. Note that throughout the 
whole paper we will use an overline to denote the time average of any operator or function. 
More precisely, if the total Hamiltonian has nondegenerate energy gaps, the average distance 
between psit) and the time averaged state is bounded by [5] 

1 ' 



where 



represents the effective dimension of the bath. Here and in the following wc characterize 
the distance between quantum states by means of the trace distance. Given two quantum 
states p^ and p^, their trace distance is defined as 

D{p\p') = l\\p'-p'l (6) 

where the trace norm is considered. The upper bound in Q, which has been proven in [5] 
for a pure product initial state, can be easily extended to a mixed product initial state, see 
Appendix A. 

On that account, the open system equilibrates under very general conditions, but 
nevertheless, despite some significant results [3l O [6l 19] , it is still a widely open question 
which are the conditions that determine whether an open system, besides equilibrating, 
thermalizes. In the following, we will focus on this issue and, in particular, on finding 
conditions on the dependence of the time averaged state on the initial state of the open 
system. 



3. Information trapping 

3.1. Time averaging map 

First of all, let us take a closer look at the time averaged state defined in In the following, 
we assume for simplicity a non degenerate Hamiltonian H — but analogous 

considerations can be done in the degenerate case. For a fully generic initial total state psBi 
one has 

p^ = Y,{Ek\psB\Ek)al, (7) 

k 
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where the notation 

4^trs{|i?fc>(^fc|} (8) 
has been introduced. If the initial total state is a product state ps fS" ps, with fixed 



environmental state pB, equation ([7|) defines a map A on the state space of the open system 
Sins): 

— f /■* 

PS Aps ■■= lim - / drps(r) V(£^fc|ps ® (9) 
that can also be written as 



t.ps^Y.Vk^s. (10) 

with 



pi.i 

k 



Pk = {Ek\ps®PB\Ek)^T:^{\Ek){Ek\ps®PB}- (11) 

This map associates to any initial state of the system the corresponding time averaged state, 
and therefore we will call it the time averaging map. This is a linear, trace preserving and 
completely positive map and its image K{S{T-Ls)) = IniA can be identified as the set of 
equilibrium states of the reduced dynamics. A state ps will be said to be invariant if it is 
left unchanged by the time averaging map, i.e. if Kps = Ps- A natural question is then 
whether equilibrium states are invariant, i.e., if, given a state ujs = ^ps for some initial 
state Ps, one has h^ujs — uis- This can happen for any ujs if and only if the map A is a 

2 2 

projector. That is, it satisfies the idempotcnce relation A = A, where of course A indicates 
the composition of A with itself. Note that 

^^Ps ^Y.^Ek\ps(^ pB\Et,){Ek'\a% ® pB\Ek')<ys- (12) 

kk' 

Let us now make the following important remark. First, recall that any trace preserving 
and positive map A is a contraction for the trace distance [TH [IS], i.e. D{Ap^,Ap^) < 
D{p^,p^) for any p^ and p^. A map is further said to be strictly contractive [ini [TT] if 
D{Ap^ , Ap'^) < D{p^,p'^) for any p^ ^ p^. Indeed, the time averaging map A is contractive, 
but in general not strictly contractive. It is clear that the only way for it to be both strictly 
contractive and idempotent is to map every initial state to the same time averaged state. 
That is, the maps Ap defined as 

Apps^P ypseSiUs) (13) 

for a fixed state p, are the only idempotent and strictly contractive maps on 5(Hs)- In fact, 
let uji,uj2 € ImA be two elements of the image of A, i.e. wi = Ap\, and u!2 = Ap^g for some 
PsjPs € '^('Hs). The idempotence of A implies that Z3(Aa;i, Aa;2) = D{uji,uj2) and, because 
of the strict contractivity, it follows that lui ^ lu^. Hence the image of A consists of only a 
single element which proves our claim. 

As a consequence, the dependence of the equilibrium state on the initial state of the 
system can be always related to the violation of either property. In other words, the absence 
of thermalization can always be associated with cither the lack of idempotence or strict 
contractivity of the time averaging map. In the following we will focus on the violation of 
the idempotence of A, which will be referred to as information trapping. This will be shown 
to capture an interesting dependence of the equilibrium state on the initial state, and we 
will demonstrate a connection between information trapping and the creation of correlations 
between the system and the bath due to their mutual interaction. 
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3.2. Measure for information trapping 

Rather than simply assessing whether the time averaging map is idempotent, one needs to 
quantify its possible deviation from idempotence in order to point out if this can be treated 
as "small". The very definition of the equilibration of a quantum system interacting with 
a finite dimensional bath involves the idea that the reduced state ps(t) will stay most of 
the time in a neighborhood of the corresponding equilibrium state Apg. For the sake of 
concreteness, let us denote as X the radius of such neighborhood. Now, if the distance 
between two different equilibrium states Ap^ and Ap| is smaller than X, the corresponding 
time evolved states Pg{t) and p'g{t) can be close to each other for almost all times, so that 
one cannot practically infer that they approach different equilibrium states by monitoring 
their evolution. This leads us to the conclusion that it is more meaningful to investigate the 
amount of information trapping of a given dynamics, rather than its mere presence. 
In particular, we propose the following measure for information trapping: 

r(A)= max D(K%sAps)- (14) 

PseS{-Hs) 

This directly quantifies the violation of the idempotence of A, and it is indeed equal to if 
and only if A is idempotent. In Appendix B, we introduce an alternative, but qualitatively 
equivalent, measure. Now, if T(A) exceeds X , there is some ps such that one can actually 
determine that ps and A/55 evolve to different equilibrium states and no thermalization 
occurs. In addition, as will be shown by means of examples, the measure T(A) provides 
a useful way to describe how the different features of a given dynamics can enhance or 
decrease the information trapping and thus the dependence of the equilibrium state on the 
initial state of the open system. 



4. Information trapping and system-environment correlations 

In this section, we explicitly connect the notion of information trapping with the interaction 
induced correlations between the system and the bath, as well as the changes in the 
environmental state. First of all, it is useful to come back to the full unitary dynamics, 
where the time averaging can be described by means of a trace preserving and completely 
positive map U, such that (compare with ([t])), 

— If* 

UpsB = lim - / drpsBir) = y^{Ek\psB\Ek)\Ek){Ek\. (15) 

t-i-oo t ^ 

Indeed, this map can be defined for any initial total state, but we will focus on the case 
PSB = Ps ® Pb, with fixed pB, to guarantee the existence of the reduced map A, which can 
be expressed as 

Aps ^tTB{U{ps® pb)] ■ (16) 

In the following diagram one can see the relation between the map U on the total system 

_ _ 2 

and both the reduced time averaging map A and its two-fold application A : 

Ps ® PB — ^SB ^U{ps® Pb) (g) PB — U [i^s ® Pb) 



tr^ 



tr 



tr, 



4 



^2 



ps > W5 = hps ojs = Aps > A PS 
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where we introduced the notation ujsb = U {ps ® Pb) to indicate the time averaged state of 
the total system, so that los = tr b^sb — ~Ps and lob — ^t^s^SB — Ab are the time averaged 
states of the system and the bath, respectively. In particular, note how the reduced map 
A is obtained after resetting the total state from ujsb to ujs ® Pb- Now, the map U o n th e 



-2 

total state is always idempotent, i.e., U = U, as can be easily checked by means of (151, 



since it amounts to a von Neumann measurement of t he to tal energy. Introducing the map 



$ = tiB o U : SCHsb) SCHs) , from the diagram (17) and the idempotence of U, one 
has A ps = ^{i^s ® Pb), while Aps = ^{ujsb)- But then, since $ is trace preserving and 
completely positive , the contractivity of the trace distance implies 

2 

D{A PS, Aps) < D{usb,^s ® Pb) < D{ujsB,i^s <^ ^^b) ^ D{pb,ujb)- (18) 

The information trapping is upper bounded by the total amount of correlations between the 
system and the bath in the total time averaged state usb plus the distinguishability between 
the time averaged state of the bath ojb and the fixed initial state pb- This means that, while 
a small amount of entanglement in the energy eigenbasis prevents a full thermalization [S], 
such phenomenon will be generally favored if the dynamics builds up a small amount of 
correlations between the system and the bath, together with small changes in the state of 
the bath. 

From a physical point of view, we can explain the connection between information 
trapping and system-environment correlations by taking advantage of the notion of 
information flow associated with the changes of the trace distance between reduced states 
in the course of time [THl [H] . The basic idea is that if there is some information trapped 
into the open system when it approaches the equilibrium, then, by resetting the system- 



environment correlations as well as the bath state to their initial condition, see (17), one 
can restart an information flow between the system and the bath, thus leading the system to 

2 

a different equilibrium state. The distinguishability between the new equilibrium state A ps 
and Aps then provides a way to quantify the information trapped into the open system due 
to the system-environment correlations and the changes in the environmental state. 

The relevance of bounds, determined by correlations in the total state as well as different 
environmental states, for the trace distance among different system states has been first 
pointed out in [20] , where the time dependence of the trace distance has been related to the 
presence of initial correlations. Here however the different system states do not correspond 
to different initial conditions, but rather to the action of distinct mappings. 

5. Examples 

5.1. Product energy eigenbasis 

As a first representative example, consider a product energy eigenbasis [3], 

H=Y. Ek,kjEk,){Ek,\® \Ek,_){Ek,\. (19) 

felfe2 

For such an Hamiltonian any reduced observable of the form A = X^fci '^fei l^fci) (£^fci I 
represents a conserved quantity. Note that a non degenerate conserved quantity on the 



open system implies a product eigenbasis of the total Hamiltonian. For i? as in (19), the 
time averaging map is not strictly contractive. In fact, one has 

UJS ^Aps ^J2^E,,\ps\E,,)\Ek,){E,,\, (20) 
fci 
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implying that if we choose as initial states two different elements of the basis {\Eki)}f^^^i ds ' 
Ps = \Ej,){Ej,\ and p| \Ei,){Ei,\, we get 

D{A\E,,}{E,,\,A\Ei,){Ei,\)=.D{\E,,}{E,,l\Ei,}{Ei,\) = 1. 

On the other hand, for a product energy eigenbasis there is no information trapping, since 
the time averaging map is a projection. Even more, as it clearly appears from the expression 



of the time averaged state (20), by setting ujs(D Pb as initial total state, the reduced system 
does not evolve at all. This clearly shows that, unlike violation of strict contractivity, 
information trapping describes a mechanism preventing thermalization which is not merely 
due to conserved quantities of the open system. 

Moreover, for a product energy eigenbasis, the total time averaged state is a product 
state, i.e. losb — ® ^b, but in general the time averaged state of the bath lob will be 
different from the initial state pB'- 

D{ub,Pb)^\\\ {EuMEkO\Ek;){Ek'^\\\. (21) 

2 

Nevertheless, we have just shown that h ps = Aps for any ps, whichever the value 



D{ujb,Pb) in equation (21). Indeed, inequality (18) gives an upper bound to the amount 
of information trapping which implies that it may well happen that, despite strong system- 
environment correlations or changes in the environmental state, the equilibrium state 
presents no information trapping. 



5.2. J aynes- Gumming s model 

Let us now consider the Jaynes-Cummings model, i.e. a two-level system interacting under 
the rotating wave approximation with a single mode of the radiation field. Moreover, the 
latter is initially in a thermal state, so that the effective dimension of the bath can be made 
arbitrarily large by properly increasing the bath temperature. Indeed, this model is much 
simpler than systems with a macroscopic number of degrees of freedom [5J |5J [7] or many- 
body systems [H [ST] |^ , which are usually taken into account when studying thermalization 
in the quantum setting. In this context, the Jaynes-Cummings model can be seen as a toy 
model, which allows us to explicitly evaluate all the quantities presented in the previous 
sections. We emphasize, however, that our general analysis can be applied to any open 
system, the only requirements being that the dimension of the total Hilbert space is finite, 
and that the open system and the bath are initially uncorrelated. 
The Hamiltonian giving the total dynamics is 

H = wo(T+(T_ + ujb^b + g{(j+®b + a^® b^) , (22) 

where (t+ = |1)(0| and (t_ = |0)(1| are the raising and lowering operators of the two-level 
system, while the creation and annihilation operators of the field mode, b^ and b, obey the 
standard bosonic commutation relation. Finally, g is the coupling constant and we will 
denote hy A ^ ujq — lo the detuning between the frequency cjg of the atom and the frequency 
uj of the field mode. Moreover, one can think of an high-energy cutoff in order to keep the 
dimension of the bath finite. For an initial total state psB — Ps ® Pb, where 

p'-iT 7) (23) 

V Poi Poo / 



Role of correlations in the thermalization of quantum systems 



8 



and pb 
by [ig 



where 



/Z is the thermal state of the bath, the reduced state at tnne t is given 



Ps{t) 



Poo{l- a{t)) + piil3{t) Piolit) 

Pml*{t) pooaW+Pii(l-/3(t)) 



(24) 



a{t) = (c'f {h,t)c{n,t))B 
I3{t) = (c^ {n + l,t)c{h + l,t))B 
-t{t) = {c{n,t)c{n+l,t))B, 
with {A) B — Tr {Apb}, the number operator n — Wh and 



c (n, t) 



COS 



v/A2 + 4g2n- - 



iA 



+ 4.g2n 

The time average can be directly calculated, thus giving 



sin ( a/A2 + Ag'^r 



POQ (1 







Poo a 



(25) 



with 



/3 = 



A2 + 4g2n 
A2 + 2g2(fj 



1) 



From (25 1 one has 



-r2 
A PS 



A2 + 4g2(n + 1) 



(A^ps)ii 



(26) 



(a'ps) 



(A'ps)ii = Poo (1 -a) (a + 13)+ pn (! + (« + /3)(/3 - 1)) 



(27) 



SO that A = A if and only ii a + /3 = 1 or a — 1 and f3 — 1. The latter case corresponds 
to 5 = 0, which implies, as expect ed f rom the discussion in section (5.1), that there is 



no strict contractivity, see equation (24 1, and A is idempotent. In all the other situations 
one has /3 + a — 1 < 1 and the map is strictly contractive. Moreover, since for A 7^ 0, 
A2+4g2fi + A2+4g2[ft^i| > 1, the only possibility to have a + /? = 1 is actually the resonant 
situation, A 



0. In this case the derivation of a through equation (26) has to be performed 
quite carefully. One can take the limit A — > into the function the series in ( 26 ) converges 
to, or equivalently note that c^(0, t)c(0, t) = 1, so that for A = one has 



1 

Z 



E 

n>0 ■ 



z 



1 

2Z 



(28) 



Thus, we have 0I + j3 ~ \ = 1/(2Z), meaning that, apart from the trivial case, in this model 
there is always information trapping and there never is only one equilibrium state. 

We can now characterize the dependence of the equilibrium state of the open system 
on its initial state by means of the measure for information trapping introduced in equation 



(14|) of section ( |3.2| ): 

r(A) 



max |pii(a + /? 

pii 

+ (1- 



-2){a + l3- 



l) + (l-a)(a + /3-l) 



(29) 
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where the maximum is assumed for ps — Note that from (28 1 one has for A = 



r(A) = ^(l-e-^") 



(30) 



In figure (|lj), we have plotted the measure T(A) as a function of the detuning A, for different 
values of the bath temperature T. The resonant situation A = represents a minimum 
for the dependence of the equilibrium state on the initial state. The residual amount of 
information trapping has to be compared with the radius of equilibration X, as discussed in 
section (3.2 1. If T(A) < X, the residual dependence of the equilibrium state on the initial 



state of the open system is not enough to recognize that ps and Aps evolve to different 
equilibrium states. On the ground of numerical simulations, we can consider the right hand 
side of ^ as an upper bound to X, see also 8J. Thus, in figure ([T]) one can see that the 
information trapping is actually larger than X for high enough detuning. 




Figure 1. Measure for information trapping T(A) defined in (291 as a function of 
A/uj for g/iu^l and I3uj = 0.003 (blue line), l3uj = 0.005 (red fine) and Puj = 0.01 
(yellow fine); the values for A = can be obtained through equation ( |30[ ), as well. 
The marks on the vertical axis give the values of the r.h.s. of Q, which upper 
bounds the value of X as discussed in the text, for A = and tor the different 
temperatures: these marks are, respectively, 0.027, 0.035 and 0.050. 



From the point of view of the Hamiltonian eigenvectors the condition A = is in fact 
very peculiar: the eigenvectors of the Jaynes-Cummings Hamiltonian, the so-called dressed 
states, reduce for A = to 

|vl'±) = -^(|0,n)±|l,n-l)), (31) 

plus the vacuum state |0,0). Every eigenvector |^^) is maximally entangled in (D^ (E) C!^, 
where (D^ is the two-dimensional subspace of Hb spanned by \n) and \n — 1). Note that, 
at resonance, both the entanglement on the energy eigenbasis and the amount of residual 



information trapping, see ( 30 1, do not depend on the coupling constant g between the system 
and the bath. Thus, for the model at hand, a high entanglement in the energy eigenbasis 
ensures an (effective) independence of the equilibrium state on the initial state of the open 
system, as one can expect from [5]. 

Now, we want to explicitly quantify the role of the system-environment correlations, as 
well as the changes in the environmental state, by means of the upper bound introduced 
in (18), i.e. D{ujsb,^s 'S^ ^b) + D{lub, Pb)- This quantity can be explicitly evaluated by 
following the same strategy employed in [24] to calculate the amount of correlations in 
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the total Gibbs state, which takes advantage of the block diagonal structure of the total 
Hamiltonian ( 22 ) with respect to the dressed states. In figure ^ one can see the measure 
T(A) as a function of the detuning A compared with the upper bound. We observe how the 
latter, despite being quite far from the actual value of the measure for information trapping, 
follows its behavior from a qualitative point of view. Indeed, as follows from the bound 



(18), small system-environment correlations and changes in the state of the bath imply a 



small amount of information trapping, and therefore an equilibrium state of the two-level 
system that hardly depends on its initial state. But for the model at hand, whenever 
strict contractivity holds, we have in addition that the more the interaction induces system- 
environment correlations and changes in the environmental state, the more the equilibrium 
state will depend on the initial reduced state ps- 



0.25 - 




Figure 2. Measure for information trapping T(A) (blue line) and D{u)sb,^s ® 
ujb) + D{ijOb,Pb) (red, dashed line) vs. A/cj for g/io — 1 and fiu) = 0.01; ujsb is 
given by U{ps ® Pb), see (171, with ps = e~^""/Z and ps ~ |1)(1| the reduced 



state maximizing the information trapping in ( 29 1 



5.3. Structured environment 

As a complementary example, we consider now the model of a system interacting with a 
structured reservoir introduced in [2^ A two-level system is coupled to two energy 
bands with the same width 5e; the energy levels in each band are equidistant and there are 
A^i {^2) levels in the lower (upper) band. The distance AE between the central levels of 
the two bands is in resonance with the free energy of the two-level system. The coupling 
constants between the two-level system and the two bands are independent and identically 
distributed complex Gaussian random variables and their overall strength is parametrized 
by a constant A. Finally, we assume that the initial state of the environment is given by a 
maximally mixed combination of the lower band levels. Thus, by means of Hilbert space 
averaging [26 or correlated projection superoperators [25 techniques, one gets the following 
equations in the Schrodinger picture for the excited state population pii(t) and the coherence 



Pii(<) = - 7Pii(*) + 7iPii(0) 
pioit) = - (zA£; + 72/2)pio(i), 



(32) 
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where 7^ = 2'KX^Ni/5e, i — 1,2, and 7 = 71 + 72- 



Ps{t) 



Pll^(7l +726 



Poie 



-(72/2-iA£;)t 



These equations are solved by 



The time averaging map corresponding to the evolution in (33 1 is, see (|9]), 



Piiz 




Poo 





Pii(l - 



(33) 



(34) 



Indeed, for 72/71 =0 this map reduces to the identity map, while in all the other situations 
it is a strictly contractive non-idempotent map. The square of the time averaging map is 
given by 



-r2 
A PS 



Pii 



7i 







Poo 





Pii(l 



7'' 



SO that the measure for information trapping defined in ( 14 ) is 



r(A) 



7i 



7? 



7 72 {Ni+N2)^ 



(35) 



(36) 



and the maximization is obtained with ps = |1)(1|- The information trapping is completely 
determined by the ratio N1/N2 and, in particular, it vanishes only in the limit N2/N1 — >■ 0, 
which corresponds to the trivial situation A = 1, or in the limit N1/N2 — > 0. 

Finally, let us present a remark about the connection between information trapping and 
the non-Markovianity of a quantum dynamics [18, .28] . Note that the relation between the 
asymptotic state of a reduced dynamics and its non-Markovianity has been studied in [29j . 
From (33 1 one can easily obtain a time-local master equation in the form 

d 



dt 



Pit) = K{t)pit) 



(37) 



to characterize the dynamics of the two-level system. The time-local generator K(t) is in 
fact simply given by |23j 

K{t) = A{t)A-\t), (38) 

so that, for the model at hand it reads 

d 



dt 



pit) = ~iAE[a+a-,pit)] 



(39) 



ri(t) 



1 



a-pit)a+ - - {a+a-,pit)} 



r2(i) [a,p(t)a, - p(i)] 



with 



ri(<)- 



727 



r2(t) = 



7172 



1 



726 T* + 71 



(40) 



71 el"* + 72 

Such coefficients are positive at every time, implying that the reduced dynamics under 
consideration is always Markovian, both in the sense that it implies a monotonic decrease 
of the trace distance in the course of time and in the sense that it is fixed by a divisible 
family of dynamical maps [T8l [28l l30l [31] . This clearly shows that one can actually have 
information trapping also in the presence of a Markovian dynamics: more generally, the 
dependence of the equilibrium state of the open system on its initial state does not provide 
a signature of non-Markovianity according to the above mentioned definitions. 
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6. Conclusions 

In this paper, we have investigated the thermahzation of finite dimensional quantum systems, 
within the framework of the theory of open quantum systems. By only assuming an initial 
product state, we have shown how one can introduce a time averaging map on the state 
space of the open system that associates to any initial state the corresponding equilibrium 
state. In this way, we could formulate relevant questions related to equilibrium properties 
of the open system in terms of suitable properties of the time averaging map. In particular, 
the dependence of the equilibrium state on the initial reduced state can be always traced 
back to the violation of at least one of the properties of strict contractivity and idempotence 
of the time averaging map. Indeed the violation of idempotence has been shown to provide 
an indication on the amount of information about the initial system state stored in the 
equilibrium state. We have therefore dubbed this violation as information trapping. It 
has been shown to be strictly connected to the interaction induced correlations between 
the system and the bath, as well as the changes in the environmental state, which keep 
track of system-environment information flow. More precisely, small system-environment 
interactions, together with small changes in the state of the bath, lead to an equilibrium 
state with a small dependence on the initial state of the open system, as quantified by means 
of the trace distance. 

Furthermore, we have introduced a measure in order to evaluate the amount of 
information trapping of a given dynamics. This provides a way to determine how the 
different features of the dynamics influence the dependence of the equilibrium state on 
the initial state of the open system and therefore how they can favor or prevent a full 
thermalization. In particular, in the Jaynes-Cummings model one can conclude that if the 
time averaging map is strictly contractive, then strong system-environment correlations and 
changes in the environmental state imply a significant dependence of the equilibrium state 
on the initial state of the open system. Indeed, it would be important to determine whether, 
or at least to what extent, this implication holds in general. 

Finally, let us note that the present results could provide a further insight into the role 
of the weak coiipliiig assumption into the process of thermalization. If the open system 
and the bath are weakly coupled, one expects that the total state at a generic time can be 
effectively described by neglecting the system-environment correlations and the changes in 
the environmental state. In this regard, it will be of interest to investigate the connection 
between the correlation properties of the total time averaged stated studied in this work and 
the correlation properties of the total state in the course of time. 
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Appendix A. Extension of the bound in the average distance between a state 
and its time average to a generic initial product state 

Here we prove the inequality Q for a generic initial product state ps ® PBi under the 
assumption that the total Hamiltonian H has non degenerate energy gaps, i.e. Ek — E^' = 
Ej — Ej> implies k — k' and j = j' or fc = j and k' = j' . We will essentially follow the proof 
in [3] for initial pure product states. First, let us introduce the notation 

Ckk' = {Ek\ps®PB\Ek'), (A.l) 

so that 

Ps{t)^J2''~'^''"~^'''^''"^k'tTB{\Ek){Ek'\} (A.2) 

kk' 

and therefore, compare with Eqs.Q and 

Aps = ^CfcfctrB{|i;fe)(£;fc|}. (A.3) 

k 

From the bound ||p|| < Vds\\p\\HS, where || • \\hs denotes the Hilbert-Schmidt norm 
II ^11 IfS ~ trsA^^, and by exploiting the concavity of the square root we have 



Dips{t),Aps) < y dstrs {(psit) ~ Aps)^}. (A.4) 



From Eqs.(A.2| and (A.3), it follows that 

ps{t)-^PS=Y. ^''^^''~''^''^'0kk'XXB{\Ek){EyW . 
k^k' 

Using the identity 

trs {trB{\Ek){Ek'\}trB{\Ek'){Ek\}}^irB{irs{\Ek){Ek\}trs{\Ek'){Ek'\}}, 
and the fact that H has non degenerate energy gaps one finds 
trs{(PsW-Ap5)2} = ^ Ckk'Ck'k^^B{^^s{\Ek){Ek\}i^s{\Ek'){Ek'\}} 

k^k' 

< J2 '"'^k' '"^''^^^B {trs {\Ek){Ek\} trs { 1^;^' ) {Ek, | }} . 

kk' 

Further exploiting the Schwarz inequality 

Ckk'Ck'k < CkkCk'k' (A. 5) 

we finally come to 



trs { {ps {t) - Aps)^} < CkkCk'k'ti-B {trs {\Ek){Ek\} trs {\Ek' > {Ek' \ }} 

kk' 

= tTBlj2^kktTs{\Ek){Ek\}J2''^k'k'tTs{\Ek'){Ek'\} \ 

{ k k' ) 

= trs [pBitf] , (A.6) 
which, together with (A.4) and ([s]), gives 
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Appendix B. An alternative measure for information trapping 



The measure for information trapping in equation ( 14 1 directly quantifies the effect of 



removing correlations and resetting the environmental state to its initial condition in the 



time averaged state ojsb, see also (17l and (18). On the other hand, it is in some sense 



arbitrary to consider only the two-fold application of the time averaging map instead of a 
high number of applications. It is in fact clear that the dynamics obtained by resetting the 
total state to Aps ® pb can still present some information trapped into the new equilibrium 
state as a consequence of further system-environment correlations built up by the interaction, 
and so on. For this reason, we introduce the following alternative measure for information 
trapping: 

roo(A)= max D( lim aVs,Aps), (B.l) 

which is set equal to 1 if the limit does not exist for some ps- It is important to note that the 
two measures, T(A) and 7^ (A), give the same qualitative characterization of information 
trapping, i.e. also 7^(A) is equal to if and only if A is idempotent. The "if part is 
obvious. To check the "only if part, 7^ (A) = implies by definition the existence of the 



limit in (B.l), and therefore in particular 



lim \\A%s ~ a' VsII =0 Vps G S{ns). (B.2) 

k—^oc 

Moreover, we can define the map A°° through A°° pg = limj._j.oo a''p5, ^^nd 7^ (A) — 
further implies that A is equal to A. This requires ImA = ImA, but then, since in 

2 k 2 k 

gener al ImA C ImA C . . . C ImA , it follows that ImA = ImA = . . . = ImA . Finally, 

T I _ _ _ _ 2 

(B.2) is equivalent to \iuik^ao A\j^-^k-i ~ 1 implying Ajj-^-^ = 1, i.e. A = A, which 
completes our proof. 

Furthermore, note that if A is strictly contractive, due to the Banach fixed point theorem 
it has a unique invariant state = limi;_j.oo A ps, which is then a natural reference state 
to quantify the dependence of the equilibrium state on the initial condition: in this case. 



the measure (B.l ) can be simply written as 



TooiA)= max D{pl,Aps). (B.3) 

psesin) 



For example, in the Jaynes-Cumniings model, for g ^ 0, the time averaging map A ( 25 ) is 
strictly contractive. Since 

r-l 

(AVs)^^= (l-a)^(^ + a-l)'+pn(^ + a-l)^ (B.4) 



1=0 



the limit map A°° is given by 



-oo 



2-/3-5 / 

which then provides the unique invariant state of the strictly contractive map A. The 




measure (B.l) is thus given by 



rc.(A)^ "+^ (B.6) 
2 — p — a 
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to be compared with (29). Analogously, in the model considered in section 5.3 for 



N2/N1 ^ 0, the map A associates every state pg with th e un ique fixed point, the vacuum 



state |0)(0|, of the strictly contractive map A in equation (34 1, i.e 



1 

and the measure for information trapping reads 



to be compared with (36). In both cases the use of this alternative measure does not quali- 



tatively change the results, which justifies to concentrate on idempotence. 
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